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Above-water  measurements  of  water-leaving  radiance  are  widely  used  for  water-quality  monitoring  and 
ocean-color  satellite  data  validation.  Reflected  skylight  in  above- water  radiometry  needs  to  be  accurately 
estimated  prior  to  derivation  of  water-leaving  radiance.  Up-to-date  methods  to  estimate  reflection  of 
diffuse  skylight  on  rough  sea  surfaces  are  based  on  radiative  transfer  simulations  and  sky  radiance  mea¬ 
surements.  But  these  methods  neglect  the  polarization  state  of  the  incident  skylight,  which  is  generally 
highly  polarized.  In  this  paper,  the  effects  of  polarization  on  the  sea  surface  reflectance  and  the  subse¬ 
quent  water-leaving  radiance  estimation  are  investigated.  We  show  that  knowledge  of  the  polarization 
field  of  the  diffuse  skylight  significantly  improves  above-water  radiometry  estimates,  in  particular  in  the 
blue  part  of  the  spectrum  where  the  reflected  skylight  is  dominant  A  newly  developed  algorithm  based 
on  radiative  transfer  simulations  including  polarization  is  described.  Its  application  to  the  standard 
Aerosol  Robotic  Network-Ocean  Color  and  hyperspectral  radiometric  measurements  of  the  1.5-year 
dataset  acquired  at  the  Long  Island  Sound  site  demonstrates  the  noticeable  importance  of  considering 
polarization  for  water-leaving  radiance  estimation.  In  particular  it  is  shown,  based  on  time  series  of  col¬ 
located  data  acquired  in  coastal  waters,  that  the  azimuth  range  of  measurements  leading  to  good-quality 
data  is  significantly  increased,  and  that  these  estimates  are  improved  by  more  than  12%  at  413  nm.  Full 
consideration  of  polarization  effects  is  expected  to  significantly  improve  the  quality  of  the  field  data 
utilized  for  satellite  data  validation  or  potential  vicarious  calibration  purposes.  ©  2012  Optical  Society 
of  America 

OCIS  codes:  010.0010,  280.0280,  010.4450,  010.1320. 


1.  Introduction 

Spectral  water-leaving  radiance,  as  measured  by 
in-water  or  above-water  instrumentation  or  from 
satellites,  is  operationally  utilized  for  ocean  and 
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coastal  water-quality  monitoring  as  well  as  in  var¬ 
ious  applications  such  as  ecological  or  climatic  mod¬ 
els,  oceanic  impacts  of  the  global  change,  fisheries, 
and  so  forth  [1].  However,  the  estimation  of  the 
water-leaving  radiance  from  satellites  remains  com¬ 
plex  due  to  the  presence  of  the  atmosphere  between 
the  water  body  and  the  sensor  [2],  As  a  result,  this 
estimation  needs  sophisticated  data  processing, 

2o \  XlXi-^VOS' 


which  is  particularly  challenging  in  coastal  water 
areas.  A  mandatory  effort  is  then  dedicated  to  vali¬ 
dation  of  the  satellite  retrievals  on  the  basis  of  com¬ 
parison  with  field-truth  measurements  [3—5]. 

The  water-leaving  radiance  can  be  derived  from 
underwater  or  above-water  instrumentation.  In  un¬ 
derwater  methods,  the  upwelling  radiance  and  down- 
welling  irradiance  are  measured  at  several  depths. 
These  data  are  then  extrapolated  to  the  surface  level 
and  propagated  upward  through  the  sea-air  inter¬ 
face.  Accuracy  in  the  retrieved  water-leaving  radiance 
depends  on  exposure  errors,  such  as  ship  shadows,  in¬ 
strument  self-shading,  biological  fouling,  and  so  forth. 
Moreover,  in  turbid  waters,  or  in  the  red  part  of  the 
spectrum,  where  the  absorption  coefficient  is  large, 
the  extrapolation  to  the  surface  may  be  particularly 
difficult  [6]. 

In  contrast,  the  above-water  method  seeks  to 
derive  the  water-leaving  radiance  from  radiometric 
measurement  of  the  sea  surface  from  the  deck  of  a 
boat  or  from  instrumentation  installed  on  a  platform. 
At  the  same  time,  the  downwelling  irradiance  is 
directly  measured  or  estimated  from  sky  measure¬ 
ments  in  order  to  normalize  the  sea  measurements 
[7,8].  The  sun  and  sky  light,  which  have  undergone 
specular  reflection  on  the  sea  surface  in  the  direction 
toward  the  radiometer,  must  be  subtracted  prior  to 
water-leaving  radiance  estimation.  This  procedure 
corresponds  to  the  most  critical  step  in  above-water 
radiometry  achievement,  as  recently  discussed  [9-11]. 
To  handle  this  issue,  it  was  first  proposed  to  remove 
the  reflected  sky  light  by  subtracting  the  sky  radiance 
measured  in  the  direction  for  which  the  sky  light 
would  be  reflected  toward  the  sensor  by  a  flat  surface 
after  multiplication  by  the  scalar  Fresnel  coefficient 
(i.e.,  first  term  of  the  Fresnel  matrix)  at  the  same  an¬ 
gle  [12].  This  approach  was  then  extended  to  account 
for  integration  of  the  scalar  Fresnel  coefficient  over 
the  radiometer  field  of  view  [13].  At  the  same  period, 
typical  values  of  the  reflection  factor  of  the  rough  sea 
surface  were  derived  from  radiative  transfer  computa¬ 
tions  [14]  and  have  been  used  to  date  in  most  of  the 
above- water  radiometiy  programs  [4,9,15].  However, 
it  was  assumed  in  all  those  methods  that  polarization 
features  occurring  during  light  reflection  on  the  sea 
surface  can  be  neglected  and  the  polarization  state 
of  the  incident  sky  light  is  not  accounted  for.  Other 
studies  made  use  of  the  polarization  properties  of 
the  specular  reflection  on  the  sea  surface  to  minimize 
the  impact  of  the  reflected  sky  light  on  the  measure¬ 
ments  by  using  a  polarizer  and  appropriate  viewing 
configurations  [16-18].  However,  no  strict  considera¬ 
tion  of  the  polarization  state  of  the  incident  sky  light 
was  taken  into  account. 

Solar  radiation  is  initially  unpolarized  when  enter¬ 
ing  the  Earth’s  atmosphere.  Yet  natural  sunlight  is 
scattered  by  aerosols  and  molecules,  refracted  and 
reflected  at  the  atmosphere-ocean  interface,  and 
further  scattered  by  hydrosols  and  water  molecules. 
As  a  result  of  these  interactions,  solar  radiation 
becomes  partially  polarized,  and  the  resulting 


polarization  field  depends  on  the  optical  properties 
of  the  atmosphere-ocean  system.  Thus,  the  downwel¬ 
ling  skylight  might  be  highly  polarized  [19],  and  its 
polarization  state  is  dependent  on  the  actual  optical 
properties  of  the  atmosphere  [20-22].  Accordingly, 
the  reflected  sky  light  on  the  sea  surface  can  be  ex¬ 
pected  to  be  itself  impacted  by  the  polarization  field 
of  the  atmosphere.  The  main  goal  of  this  study  is 
then  to  investigate  the  effects  of  polarization  on  the 
sea  surface  signal  and,  in  turn,  on  the  performances 
of  the  water-leaving  radiance  estimation  from  above¬ 
water  measurements. 

This  paper  is  organized  as  follows.  The  following 
section  describes  the  radiative  transfer  modeling  of 
the  sky  light  reflection  on  a  rough  sea  surface  and 
theoretically  discusses  the  impact  of  polarization  on 
the  estimation  of  the  surface  reflectance  and  water¬ 
leaving  radiance.  The  third  section  summarizes  the 
characteristics  of  the  Long  Island  Sound  Coastal 
Observatory  (LISCO),  which  combines  collocated 
multi-  and  hyperspectral  radiometers  installed  on  a 
platform.  Measurements  of  the  sky  and  sea  radiances 
and  of  the  downwelling  irradiance  have  been  acquired 
since  October  2009  for  water-leaving  radiance  estima¬ 
tion  and  aerosol  determination  via  the  Aerosol 
Robotic  Network  (AERONET)  program  [23].  The 
fourth  section  describes  the  polarization-based  correc¬ 
tion  scheme  for  water-leaving  radiance  retrieval  from 
above-water  measurements.  Results  of  applications  to 
the  LISCO  data  are  presented  and  discussed  based  on 
comparison  of  collocated  field  data.  Conclusive  re¬ 
marks  and  perspectives  are  given  in  the  last  section. 

2.  Theoretical  Background 


A.  Above-Water  Signal 

The  total  bidirectional  reflectance  just  above  the  sea 
surface  is  defined  at  a  given  wavelength  as  follows: 


PtiPs<0vi4>) 


xLt{0s,0y,<l>) 

'EM) 


(1) 


where  Lt  is  the  upward  radiance  measured  from  an 
above-water  radiometer  pointing  at  the  sea  with  a 
given  viewing  configuration  defined  by  the  solar 
zenith  angle  08 ,  the  viewing  zenith  angle  0Vf  and  the 
relative  azimuth  between  the  Sun  and  the  radio¬ 
meter  cj).  E8  is  the  downwelling  total  irradiance  at 
the  sea  level  just  above  the  surface.  In  like  manner, 
the  definition  of  Eq.  (1)  can  be  extended  to  the  Stokes 
vector,  expressed  in  reflectance  unit,  S  =  \p,q,  u,v]T, 
which  describes  electromagnetic  radiation,  including 
its  polarization  state,  in  terms  of  directly  measurable 
quantities;  the  superscript  T  stands  for  the  trans¬ 
pose  of  the  vector. 

The  Stokes  vector  is  defined  in  this  paper  with 
respect  to  the  reference  frame  determined  by  the 
viewing  direction  of  the  above-water  radiometer  and 
the  local  zenith  (i.e.,  the  local  meridian  plane).  For 
incoherent  light,  Stokes  vectors  are  additive,  which 
enables  us  to  decompose  the  above-water  signal  at 
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a  given  wavelength  and  a  given  viewing  configura¬ 
tion  as  follows: 


St  =  4-  S8urf  4*  4*  Sy,.  (2) 

Sf ,  Sgt  S8Urf ,  S^,  and  S^,  stand  for  the  Stokes  vectors  of 
the  total,  sun  glint,  sea  surface  reflection  of  sky  light, 
foam,  and  water-leaving  components  of  the  above¬ 
water  signal,  respectively.  More  precisely,  is  the 
contribution  arising  from  specular  reflection  of  direct 
sunlight  on  the  sea  surface,  Ssurf  is  the  contribution 
generated  by  specular  reflection  of  atmospherically 
scattered  light  on  the  sea  surface,  is  the  contribu¬ 
tion  arising  from  sunlight  and  sky  light  reflected  by 
whitecap-covered  areas  of  the  sea  surface,  and  Sw  is 
the  contribution  of  underwater  light  scattered  to¬ 
ward  the  above-water  radiometer  after  propagation 
through  the  sea  surface. 

Generally,  the  data  processing  of  above-water 
measurement  permits  filtering  out  any  data  con¬ 
taminated  by  prohibitive  amounts  of  sun  glint  or 
foam  contribution  [4,24] .  Consequently,  this  study 
will  be  mostly  dedicated  to  analyzing  the  surface 
reflection  of  diffuse  skylight  and  the  water-leaving 
contribution,  and  Eq.  (2)  can  be  reformulated  as 


S,  =  Ssurf  4-  4-  e ,  (3) 

where  the  term  €  accounts  for  residual  sun  glint  or 
foam  that  might  remain  even  after  data  filtering.  In 
neglecting  e ,  the  main  issue  to  solve  in  above-water 
radiometry  is  to  accurately  estimate  S8Urf  in  order  to 
retrieve  from  measurements. 

B.  Fresnel’s  Reflection  Matrix 

The  Stokes  vector  of  the  surface  contribution  just 
above  water  can  be  first  calculated  for  a  flat  surface. 
This  computation  is  based  on  Fresnel's  reflection 
matrix,  for  which  a  beam  of  light  is  reflected  toward 
the  sensor.  Let  n  be  the  complex  refractive  index  of 
sea  water  relative  to  the  air  index  and  0,*  be  the  in¬ 
cident  angle  of  a  beam  of  light  illuminating  the  sur¬ 
face,  which  is  equal  to  half  the  phase  angle  (i.e,, 
complementary  of  the  scattering  angle).  The  angle 
0,  can  be  defined  with  respect  to  the  above-water 
viewing  configuration  for  any  orientation  of  the  flat 
surface  considered: 

cos  2 0,  =  -  cos  O'  cos  6V  -  sin  O'  sin  0V  cos  <f>' ,  (4) 

with  ff  and  <p'  the  zenith  angle  and  the  relative 
azimuth  with  the  viewing  direction  of  the  incident 
beam  considered,  and  0V  the  viewing  zenith  angle. 
For  more  details  on  the  geometrical  parameters 
used,  the  reader  is  referred  to  Fig.  1.  Note  that  the 
relative  azimuth  <p'  is  equal  to  0°  when  the  source 
of  the  light  beam  being  considered  and  the  sensor 
are  in  opposition,  following  the  convention  commonly 
used  in  above-water  radiometry. 

For  a  given  incident  angle,  the  Fresnel's  reflection 
matrix  is  expressed  in  the  plane  defined  by  the  Sun 


Fig.  1.  Definition  of  the  viewing  geometry  convention  used:  0i  is 
the  incident  angle  of  a  beam  of  light  illuminating  the  surface, 
0  and  <f>  are  the  zenith  angle  and  the  relativo  azimuth  with  the 
viewing  direction  of  the  incident  beam  considered,  0V  is  the  viewing 
zenith  angle,  and  0  is  the  scattering  angle. 


and  satellite  directions,  that  is,  the  plane  of  inci¬ 
dence,  as  follows: 
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where  the  superscript  *  stands  for  the  complex  con¬ 
jugate.  The  terms  rx  and  rB  are  the  components  per¬ 
pendicular  and  parallel  to  the  plane  of  incidence, 
respectively,  defined  by 


_  cos  0,  -  y/n2  -sin20, 

cos  Oj  +  y/n2  -  sin2  0,- 
_  n2  cos  Oi  -  y/n2  -  sin2  Oi 
11  n2  cos  0{  +  y/n2  -  sin2  0, 


(6) 


The  ry  matrix  has  been  defined  with  respect  to  the 
incident  plane.  In  order  to  formulate  the  Stokes 
vector  originating  from  specular  reflection  with  re¬ 
spect  to  the  reference  plane,  the  tf  matrix  must  be 
reprojected  accordingly  to  get  the  proper  reflection 
matrix  fy: 


=  L(*  ~  (7) 


where  L  is  the  rotation  matrix  for  the  angles  a\  and 
a2  defined  with  respect  to  the  solar  and  viewing 
zenith  angles  and  the  scattering  angle  0  (given  by 
n-2  0,): 
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COS  (5\  = 

cos  <r2  = 


cos  6V  -  cos  0  cos  0 
sin  0  sin  0 
cos  0  -  cos  0V  cos  0 
sin  0V  sin  0 


(8) 


In  the  case  of  a  flat  surface,  the  Stokes  vector  Sg^f 
can  be  therefore  expressed  as 

Ssurf 4>)  -  ^W)Ssky(^, 0  =  x-Ou,4>),  (9) 

where  Ssky  is  the  Stokes  vector  of  the  diffuse  sky  light 
as  measured  by  a  sensor  looking  upward  with  an  an¬ 
gle  0V.  As  a  result,  the  component psurf ,  corresponding 
to  the  reflectance  of  the  sea  surface,  is  obtained  by 

Amrf  =  R\\P*Vy  +  ^12<7sky  +  ^13wsky  +  ^H^sky*  (10) 

where  the  terms  Rn  to  R14  are  the  components  of 
the  first  row  of  the  Fresnel  matrix.  It  is  worth  noting 
that  the  fourth  component  t>gky  of  Ssky  is  very  low  and 
can  generally  be  neglected  in  comparison  to  the  first 
three  components  [25-27].  The  respective  values  of 
i?ii,  Ri2 ,  and  Rx 3  are  therefore  of  paramount  impor¬ 
tance  in  the  generation  of  the  upwelling  signal  of  the 
sea  surface  contribution  as  measured  by  above-water 
instrumentation. 

The  Fresnel  matrix  component  Rn  and  the  compo¬ 
nents/?^  and/?i3  normalized  by  Rn  are  displayed  in 
Fig.  2  for  a  given  viewing  angle  equal  to  40°  and  solar 
angle  ranging  from  0°  to  70°  for  the  full  azimuth 
range.  It  can  be  readily  seen  in  this  figure  that  the 
term  Rn  is  virtually  constant  and  equal  to  0.021  for 
a  large  range  of  azimuth  roughly  comprised  between 
-60°  and  60°.  Regarding  the  rest  of  the  azimuth 
range,  this  term  is  strongly  dependent  on  the  solar 
angle  and  can  reach  values  up  to  0.04.  More  interest¬ 
ingly,  the  contribution  of  /?12  and  Rx^  are  far  from 
being  negligible  in  comparison  to  Rn*  In  particular, 
the  values  of  these  two  terms  can  correspond  to  up  to 
50%  of  the  i?n  values,  as  shown  by  the  ratios  Rn/Rn 
and  Ris/Rn,  for  almost  the  whole  range  of  azimuths. 
In  conclusion,  the  nondiagonal  terms  of  the  Fresnel 


matrix  must  be  taken  into  account  to  accurately  com¬ 
pute  the  bidirectional  reflectance  of  the  sea  surface, 
even  for  a  flat  surface. 

In  general,  the  sea  surface  is  roughened  by  the 
wind.  As  a  result,  a  myriad  of  reflecting  facets  might 
reflect  skylight  originating  from  any  direction  of  the 
sky.  Various  contributions  of  the  skylight  will  reach 
the  above-water  sensor  due  to  specular  reflection 
from  waves  and  wavelets  at  the  sea  surface.  Follow¬ 
ing  Cox  and  Munk  [28-30],  the  ocean  surface  can  be 
modeled  based  on  a  distribution  of  small  facets  that 
are  oriented  following  a  near-Gaussian  distribution. 
In  this  model,  the  statistical  distribution  of  the  orien¬ 
tation  of  the  facets  leading  to  sunlight  reflection  to¬ 
ward  the  radiometer  is  given  with  respect  to  the 
viewing  configuration  as  well  as  the  wind  speed  and 
azimuth,  denoted  as  ws  and  <£wind>  respectively.  Thus, 
the  Stokes  vector  of  the  sea  surface  contribution  S8urf 
can  be  written  as  (see  [31],  for  instance) 


S*urf  = 


Rf(0j) 

4coaOvcosAQn 


p{0,Ov.4>f\w8,4>mirA) 


xSBky(0.,^dW. 


(11) 


The  term  0n  is  the  angle  between  the  zenith  and  the 
normal  to  the  facets,  and  is  computed  from  Eq.  (4). 
The  term  p  is  the  wave  slopes  distribution,  and  C, 
accounts  for  the  shadowing  effects  of  the  waves  [32] 
and  depends  on  the  viewing  geometry  and  the  wind 
speed  (see  [31]  for  the  complete  expression  of  this 
term).  To  solve  Eq.  (11)  a  complete  knowledge  of  the 
term  Ssky,  for  any  incident  direction,  is  needed. 
Recent  instrumental  developments  have  enabled 
us  to  measure  this  term  for  the  radiance  parameter 
(i.e.,  first  component  of  the  Stokes  vector)  and  the 
polarization  components  [19,33],  but  further  efforts 
are  still  to  be  done  for  their  operational  use  in 
above-water  radiometry.  A  practical  solution  to  fully 
derive  Ssurf  is  therefore  to  solve  the  radiative  transfer 
equation,  which  can  be  easily  achieved  based  on  the 
numerous  codes  available,  including  polarization 
computation;  see  [27,31,34-37],  to  cite  only  a  few. 


R,t  Rl2fRn  Ru/Rit 

0.02  0.03  0.04  - 1.0  - 0.5  0,0  0.5  -1.0  -OS  0.0  OS  1.0 


O'  O'  O' 

Fig.  2.  (Color  online)  Fresnel  coefficient  Rn  and  the  Fresnel  coefficients  R12  andi?l3  normalized  byRu.  The  viewing  zenith  angle  is  set  to 
40®  The  concentric  circles  represent  the  solar  zenith  angles  by  step  of  10°  (from  0°  to  70°),  and  angles  in  the  polar  diagrams  represent  the 
relative  azimuth  with  the  Sun.  By  convention,  the  azimuth  is  equal  to  0°  when  the  Sun  and  the  sensor  are  in  opposition. 
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C.  Estimation  of  the  Water-leaving  Reflectance 

In  above-water  radiometry,  the  standard  way  to  de¬ 
rive  the  sea  surface  reflectance  is  based  on  sky  radi¬ 
ance  measurements  Lsk y  acquired  at  the  same  time 
as  the  total  sea  radiance  Lt  [4,7,38].  In  order  to  re¬ 
move,  or  at  least  minimize,  the  impact  of  the  sun 
glint  on  the  measured  sea  radiance,  the  average  L J 
is  computed  over  a  fixed  percentage,  typically  5%,  of 
the  lowest  values  of  sea  radiance  measured  during 
the  acquisition  sequence  [7,38].  The  water-leaving 
radiance  Lw  is  then  computed  as 

Lw (0$ <0V %<f>)  as  L; ( 08 ,0v,<f>)-f(0s,0v,<p,  ws) 

x  L8ty(08,  ff  =  n  -  0Vi  <p).  (12) 

The  second  term  of  the  right-hand  side  of  the 
equation  is  the  estimate  of  the  radiance  resulting 
from  diffuse  skylight  reflected  off  the  sea  surface. 
The  coefficient  f(0s,0v,<p,ws)  is  the  sea  surface- 
reflectance  factor,  which  is  dependent  on  sky  condi¬ 
tion,  wind  speed  ws ,  and  the  viewing  configuration. 
This  term  was  parameterized  from  radiative  transfer 
computation  in  neglecting  the  incident  polarization 
field  [14].  The  water-leaving  reflectance,  pw ,  is  then 
obtained  by  multiplying  Lw  by  n/Er  In  this  study,  we 
propose  to  replace  the  estimate  of  the  sea  surface 
signal  by  that  computed  through  Eq.  (11),  which 
includes  the  impact  of  the  polarization  of  sky  light 
on  the  estimation  of  the  water-leaving  signal. 
Equation  (12)  can  be  then  rewritten  as 

pw(0..  ou,4>)=  *L'(6‘Ey  H  eu,4>)-  (13) 

The  sea  surface  reflectance  psurf  was  simulated 
using  Chowdhary’s  radiative  transfer  code  [341  for 
a  viewing  angle  of  40°,  a  solar  zenith  angle  ranging 
from  0°  to  90°  by  a  step  of  about  1°  for  any  azimuth 
value.  In  these  computations,  the  wave  slope  distri¬ 
bution  based  on  Cox  and  Munk  parameterization 
was  assumed  including  modeling  of  wave  self¬ 
shadowing  [32].  In  addition,  these  simulations  were 
achieved  following  two  modalities:  (1)  in  fully  consid¬ 
ering  the  impact  of  the  polarization  field  occurring  at 
the  sea  surface,  referred  to  as  p8Urf(vecf°r)>  and  (2) 
in  neglecting  the  polarization  components  of  the 
Fresnel  matrix  (i.e.,  i?12  =  i?i3  =  0),  referred  to  as 
Pgurf(scalar).  Note  that  under  this  assumption,  the 
term  pgurf  (scalar)  is  equivalent  to  the  term  /X8ky  of 
Eq.  (12)  used  in  the  current  correction  methods.  In 
comparing  p8 urf  (vector)  with  pBUrf  (scalar),  the  impor¬ 
tance  of  considering  polarization  in  the  estimation  of 
the  water-leaving  reflectance  can  thus  be  inferred. 

The  simulations  of  pBUrf(scalar)  and  p8Urf(vector) 
at  550  nm  are  shown  in  Figs.  3(a)  and  3(b)  for  two 
different  azimuths,  namely  90°  and  135°,  and  three 
different  aerosol  optical  thicknesses:  0,  0.1,  and  0.2. 
Note  that  the  microphysical  properties  of  the  aero¬ 
sols  also  influence  the  signal  reflected  by  the  sea 
surface.  For  the  sake  of  simplicity  the  results  are 
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Fig.  3.  (Color  online)  Sea  surface  reflectance  at  550  nm  for 
no-wind  conditions  with  respect  to  the  solar  angle  computed 
(a)  in  ignoring  the  polarization  terms  of  the  Fresnel  matrix  (i.e., 
R12  =  /?i3  =  0)  and  (b)  in  fully  considering  polarization  behavior 
of  skylight  reflection  at  the  sea  surface,  (c)  Relative  difference  ex¬ 
pressed  in  percent  between  the  sea  surface  reflectances  displayed 
in  (a)  and  (b).  The  results  are  shown  for  two  different  azimuths, 
namely  90°  (solid  curves)  and  135°  (dashed  curves),  and  for  differ¬ 
ent  aerosol  optical  thicknesses:  0  (blue),  0.1  (green),  and  0.2  (red). 

uniquely  shown  for  a  fine-mode  aerosol  model  char¬ 
acterized  by  a  mean  radius  of  0.17  pm  and  a  refrac¬ 
tive  index  of  1.45.  However,  similar  conclusions  could 
be  drawn  irrespective  of  the  aerosol  model  consid¬ 
ered.  When  no  polarization  is  considered  in  the 
computations  [Fig.  3(a)],  pBlirf (scalar)  is  virtually  con¬ 
stant  with  respect  to  the  solar  angle  up  to  an  angle  of 
around  50°  when  an  azimuth  of  90°  is  considered. 
Beyond  this  solar  angle,  the  estimated  contribution 


8328  APPLIED  OPTICS  /  Vol.  51,  No.  35/10  December  2012 


of  the  sea  surface  increases.  In  contrast,  when  polar¬ 
ization  is  considered  [Fig.  3(b)],  Amrf (vector)  starts 
decreasing  from  solar  angles  around  15°  and  keeps 
on  decreasing  at  important  solar  angles.  The  beha¬ 
vior  of  the  sea  surface  reflectance  with  respect  to 
the  solar  angle  is  therefore  strongly  impacted  by  the 
polarization  state  of  the  sky  light  beam  hitting  the 
sea  surface. 

A  more  quantitative  view  of  the  surface  polariza¬ 
tion  effects  can  be  obtained  by  computing  the  relative 
difference  between  /Wf  (scalar)  and  Amrf (vector)  as 
follows: 

W  =  tMft-.fr*1") I -^><vector)% 
/W(vector) 

This  relative  difference  is  shown  in  Fig.  3(c).  For  an 
azimuth  of  90°,  which  generally  is  used  in  above¬ 
water  radiometry  [4],  the  exhibited  discrepancies 
between  p8urf  (scalar)  and  (vector)  range  between 

-20%  at  smaller  solar  angles  up  to  over  +50% 
beyond  solar  angles  of  60°.  This  is  consistent  with 
expectations  as  the  polarization  terms  i?12  and  i?13 
of  the  Fresnel  matrix  can  attain  50%  of  the  Rn  value, 
as  discussed  before  [see  Fig.  (2)].  Moreover,  it  is 
worth  noticing  that  the  downwelling  sky  radiance  is 
in  general  highly  polarized,  with  a  degree  of  polari¬ 
zation  that  can  exceed  60%  for  a  large  range  of  geo¬ 
metries  in  a  purely  molecular  atmosphere  [39,40].  As 
a  result,  even  though  aerosols’  presence  tends  to 
diminish  the  degree  of  polarization,  the  high  degree 
of  polarization  of  the  downwelling  sky  light  enhances 
the  importance  of  considering  the  full  Fresnel  matrix 
for  estimation  of  the  sea  surface  reflectance. 

In  contrast,  the  observed  discrepancies  between 
/?8Urf(scalar)  and  p8Urf(vector)  are  significantly  smal¬ 
ler  when  measurements  are  acquired  at  an  azimuth 
of  135°.  For  this  viewing  configuration,  the  relative 
difference  A p8Urf  remains  negative  and  comprises 
between  -25%  and  0%  over  the  full  range  of  solar 
angles.  On  the  other  hand,  it  can  be  readily  seen 


in  Fig.  3(c)  that  the  effect  of  the  aerosol  loading  is 
largely  smoother  than  that  of  the  azimuth  used  for 
observation.  However,  the  aerosol  optical  thickness 
induces  significant  variation  in  the  absolute  surface 
reflected  signal  [see  Fig.  3(b)]  and  must  also  be  taken 
into  account  in  the  derivation  of  the  sea  surface  re¬ 
flectance.  In  other  works,  the  errors  originating  from 
neglecting  polarization  were  shown  at  the  top-of- 
atmosphere  level  [31,41].  The  results  of  these  studies 
showed  similar  dependencies  of  the  errors  on  the 
viewing  geometry  as  in  Fig.  3(c),  which  corroborates 
that  those  errors  stem  from  the  polarization  effects 
occurring  at  the  air-sea  interface. 

A  synoptic  view  of  the  impact  of  atmospheric  polar¬ 
ization  on  the  sea  surface  reflectance  can  be  obtained 
by  plotting  the  relative  difference  A />surf  for  all  the 
azimuth  and  solar  angles,  and  for  the  viewing  angle 
of  40°  commonly  used  in  above-water  radiometry;  see 
Fig.  4.  In  Fig.  4(a),  calculations  were  done  at  550  nm 
for  an  aerosol  optical  thickness  of  0.1  and  for  flat  sea 
surface.  The  results  displayed  in  this  figure  corre¬ 
spond  to  the  case  previously  shown  in  Fig.  3(c)  but 
extended  to  all  azimuth  conditions.  For  instance,  it 
can  be  seen  in  this  figure  that  the  relative  difference 
at  90°  azimuth  is  first  negative  for  solar  angles  smal¬ 
ler  than  35°  and  then  becomes  positive  for  greater 
solar  angles.  It  can  also  be  noted  that  the  large  ne¬ 
gative  differences  occur  in  the  vicinity  of  40°  solar 
angle  and  0°  azimuth.  However,  this  region  corres¬ 
ponds  to  the  zone  of  influence  of  the  sun  glint  and 
cannot  be  used  for  above-water  radiometry  purposes. 

For  wind  speeds  of  4  and  10  m  s'1  [Figs.  4(b) 
and  4(c)],  the  sea  surface  ruffles  and  becomes  a  rough 
surface  on  which  myriads  of  wave  facets  can  reflect 
light  coming  from  various  directions  of  the  sky  re¬ 
gion,  thereby  creating  depolarization  of  the  signal  by 
beam  superposition.  As  a  result,  it  could  be  expected 
that  the  impact  of  sky-light  polarization  on  surface 
reflection  is  smoothed  away  when  wind  speed  in¬ 
creases.  Interestingly,  this  behavior  does  occur,  but 
the  relative  difference  Apsurf  remains  significant 
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Fig.  4.  (Color  online)  Relative  difference  A /?<urf  at  550  nm  and  an  aerosol  optical  thickness  of  0.1  for  different  wind  speed  conditions:  (a)  0, 
(b)  4,  and  (c)  10  ms-1.  The  polar  diagrams  account  for  all  the  azimuth  and  solar  angles  for  a  given  viewing  angle  of  40p.  The  concentric 
circles  represent  the  solar  zenith  angles  by  a  step  of  10°  (from  0°  to  70°),  and  angles  in  the  polar  diagrams  represent  the  relative  azimuth 
with  the  Sun. 
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over  the  viewing  configurations.  Based  on  this  theo¬ 
retical  analysis,  it  can  therefore  be  safely  concluded 
that  the  polarization  state  of  the  sky  light  reflected 
on  the  sea  surface  must  be  accurately  taken  into  ac¬ 
count  for  estimating  the  sea  surface  reflectance  and 
in  turn  the  water-leaving  radiance. 

3.  LISCO 

A.  LISCO  Site 

The  LISCO  platform  is  located  at  around  3  km  from 
the  shore  of  Long  Island  near  Northport,  New  York, 
USA.  The  coordinates  of  the  site  are  N  40°57'16", 
W  73°20'30'\  The  bathymetry  in  the  immediate  vici¬ 
nity  of  the  platform  exhibits  a  plateau  at  around 
13  m  depth.  LISCO  is  located  in  the  area  of  western 
Long  Island  Sound,  which  is  usually  moderately 
turbid,  very  productive,  and  principally  driven  by 
phytoplankton  biomass  and  associated  detrital  ma¬ 
terials  rather  than  suspended  sediments  [42].  Steady 
surface  currents  in  the  LISCO  area  are  around 
0.3  ms"1  in  average  according  to  National  Oceanic 
and  Atmospheric  Administration  radar  data  sets. 
The  instruments  are  positioned  on  a  retractable  tower 
on  the  platform  with  an  elevation  of  12  m. 

B.  Multispectral  and  Hyperspectral  LISCO  Data 

The  platform  combines  a  multispectral  SeaPRISM 
system  (CIMEL  Electronique,  France)  with  a  collo¬ 
cated  hyperspectral  HyperSAS  system  (Satlantic, 
Canada).  The  SeaPRISM  system  is  part  of  AERONET 
and  particularly  of  its  ocean  color  (OC)  component 
[4,8,23].  The  instruments  are  positioned  on  a  retract¬ 
able  tower  on  the  LISCO  platform  with  an  elevation  of 
12  m.  Both  instruments  were  installed  in  October 
2009  and  have  been  providing  data  since  then. 

The  SeaPRISM  system  is  made  up  of  CE-318  sun 
photometers  modified  to  meet  requirements  for 
above-water  radiometry.  The  photometers  perform 
radiance  measurements  of  the  total  radiance  of  the 
sea,  Lt(08. 0v.<p ),  and  the  sky,  Ls^y(0s,  ff%  </))>  for  the  re¬ 
lative  azimuth  angle  with  respect  to  the  sun  (p  and 
the  respective  viewing  angles  0V  and  & ,  with  0V  = 
n  -  0.  Thanks  to  the  rotating  feature  of  SeaPRISM, 
the  azimuth  <p  is  fixed  to  90°  regardless  of  the  sun 
position;  the  downward  zenith  viewing  angle  6V  is 
set  to  40°.  The  spectral  bands  of  the  SeaPRISM 
system  are  centered  on  413,  442,  491,  551,  668,  870, 
and  1018  nm. 

The  collocated  hyperspectral  measurements  are 
made  by  a  HyperSAS  system,  providing  high- 
precision  hyperspectral  above-water  measurements 
of  downwelling  and  upwelling  spectral  radiances  as 
well  as  downwelling  spectral  irradiance.  The  radiance 
and  irradiance  measurements  of  HyperSAS  are 
carried  out  for  180  spectral  channels  regularly  spaced 
between  305  and  905  nm.  HyperSAS  has  two  radiance 
sensors,  one  looking  down  at  the  water  and  the  other 
looking  skywards.  These  sensors  provide  the  sea 
Lt(A,0,  <p)  and  the  sky  radiance  L8^y(X,  0<(p)  for  the 
same  viewing  angle  as  SeaPRISM,  that  is,  0V  =  40°. 


However,  the  sensors  are  pointing  exactly  westwards, 
resulting  in  a  relative  azimuth  <f>  changing  with 
respect  to  the  sun  position.  Each  HyperSAS  and 
SeaPRISM  sea-viewing  measurement  sequence  is 
executed  every  30  min  within  plus  or  minus  4  h  of 
12:00  local  time  (LT). 

The  calibration  of  the  SeaPRISM  sun  photometer 
was  carried  out  by  the  NASA  AERONET  group 
in  accordance  with  the  standard  procedures  of 
AERONET-OC.  The  HyperSAS  system  calibration 
was  carried  out  by  Satlantic,  Inc.  (Halifax,  Canada), 
and  checked  at  the  City  College  of  New  York  (CCNY) 
Optical  Remote  Sensing  Lab.  The  recalibration  by 
Satlantic,  Inc.,  of  the  HyperSAS  system  has  shown 
a  radiometric  stability,  over  a  time  period  of  approxi¬ 
mately  15  months,  of  better  than  1%  for  up-  and 
downwelling  radiance  sensors  and  better  than  0.5% 
for  the  irradiance  sensor.  In  addition,  errors  due  to 
nonperfect  cosine  response  of  the  irradiance  sensor 
were  estimated  to  be  smaller  than  2%  on  average. 
SeaPRISM  data  are  transferred  by  a  satellite  link 
to  NASA,  processed  by  the  NASA  AERONET  group, 
and  posted  on  the  NASA  AERONET  website.  The 
near-real-time  transmission  of  HyperSAS  data  is 
achieved  by  broadband  cellular  service  to  the  CCNY 
server. 

In  this  study,  the  SeaPRISM  data  used  corre¬ 
sponds  to  the  Level  2.0  data  archive.  Those  data  are 
cloud  screened  and  quality  assured,  with  up-to-date 
calibration.  The  HyperSAS  data  used  corresponds  to 
data  which  have  successfully  passed  the  data-quality 
check  as  defined  in  [9]. 

C.  Aerosol  Retrievals 

In  addition  to  these  OC  measurements,  the  regular 
data  acquisitions  of  AERONET  are  also  carried 
out,  which  permits  accurate  retrievals  of  the  aerosol 
optical  thickness  and  the  fine-coarse  aerosol  mode 
fraction  [43,44].  The  aerosol  component  of  AERONET 
[23,45]  is  a  well-established  network  of  over  200  semi¬ 
permanent  ground-based  sunphotometers  world¬ 
wide.  They  provide  standardized  high  quality  aerosol 
measurements  which  are  used  here  as  reference 
data.  Visible  and  near-infrared  measurements  of  the 
aerosol  optical  depth  are  available  with  an  accuracy 
of  about  0.01  to  0.02.  The  LISCO  Level  2.0  data  ar¬ 
chive  of  aerosol  inversion  is  used  in  this  study  with, 
specifically,  the  respective  spectral  aerosol  optical 
thicknesses  of  the  fine-  and  coarse-mode  aerosols  as 
derived  by  the  spectral  deconvolution  method  of  [44]. 

4.  Retrieval  of  Water-Leaving  Radiance  from  Above- 
Water  Measurements 

A.  Modeling  of  the  Sky  Glint  Contribution 

The  aerosols  are  assumed  to  follow  a  bimodal  size 
distribution.  As  pointed  out  by  Whitby  [46]  and 
Junge  [47,48],  a  realistic  size  distribution  can  be 
modeled  based  on  the  sum  of  single  modes.  Several 
studies  showed  that  a  satisfactory  simulation  of  aero¬ 
sols  can  be  obtained  using  a  mixture  of  a  fine  mode 
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with  a  coarse  mode  [49-53].  The  size  distribution 
N(r)  of  these  two  modes  can  be  parameterized  with 
the  same  log-normal  expression,  as  follows: 

rfN(r)  =  77T* exp  HCV-8^ loe  r' 

(15) 

where  r,  the  mean  radius,  and  a,  the  standard  devia¬ 
tion  of  log  r,  are  specific  of  a  given  aerosol  mode. 
Based  on  this  log-normal  size  distribution,  the  opti¬ 
cal  properties  of  these  two  aerosol  models  were  com¬ 
puted  through  Mie  theory  calculation  assuming 
spherical  particles.  The  microphysical  properties 
were  chosen  to  be  as  close  as  possible  to  the 
AERONET  inversion  data  typically  retrieved  at  the 
LISCO  site.  The  fine  mode  corresponds  to  a  mean 
radius  of  0.17  pm,  a  refractive  index  of  1.45,  and  a 
standard  deviation  of  the  size  distribution  of  0.46. 
The  coarse  mode  corresponds  to  a  mean  radius  of 
2.0  pm,  a  refractive  index  of  1.35,  and  a  standard  de¬ 
viation  of  the  size  distribution  of  0.75.  The  fine-mode 
aerosol  model  is  assumed  to  be  slightly  absorbing  ac¬ 
cordingly  to  the  AERONET  retrievals  [11]  with  an 
imaginary  part  of  the  refractive  index  equal  to  0.005. 

Based  on  the  respective  microphysical  properties 
of  the  fine-  and  coarse-mode  aerosols,  the  sea  surface 
reflectance  pBnr f  was  computed  through  radiative 
transfer  computation  including  polarization  effects 
as  discussed  in  Section  2.  This  reflectance  was  com¬ 
puted  for  the  fine  mode  and  coarse  mode,  separately, 
for  aerosol  optical  thicknesses  ta  ranging  from  0  to 
0.4  by  a  step  of  0.05.  The  radiative  transfer  computa¬ 
tions  were  carried  out  for  viewing  and  solar  zenith 
angles  ranging  from  2°  to  90°  with  increments  of 
around  1.5°.  The  sea  surface  reflectance  is  exactly  re¬ 
composed  for  any  azimuth  angle  based  on  a  Fourier 
series  approach.  Then,  the  total  sea  surface  reflec¬ 
tance  of  the  fine  and  coarse  aerosol  mixture  is  com¬ 
puted,  as  shown  by  Wang  and  Gordon  [54],  for  any 
viewing  configuration  defined  by  (08,0V%  <p): 

P8xxrf(TaW-.0C'0v.<t>)  =  YpfS{lr((Ta(X);0,,  0V,  <f>) 

+  (1-rKurf  (*«W);0,.<W). 

(16) 

where  (fs  .  and  pP  rf  stand  for  the  sea  surface  reflec- 
tance  calculated  for  the  fine  mode  and  the  coarse 
mode,  respectively.  They  are  calculated  for  any  aero¬ 
sol  optical  thickness  by  linear  interpolation  within 
the  exact  simulation  generated  for  discrete  xa  values. 
The  term  y  is  the  mixing  coefficient  between  the  two 
aerosol  modes  (i.e.,  fine  and  coarse  mode)  given  by 


It  is  worth  mentioning  that  Eq.  (16)  fails  to  accurately 
represent  the  Stokes  vector  of  aerosol  mixtures 
when  strongly  absorbing  aerosols  are  present. 


However,  such  an  equation  is  frequently  used  in 
OC  and  aerosol  remote  sensing  [55,56]  because  it 
permits  introduction  of  a  continuous  variable,  namely 
the  mixing  coefficient  y ,  thereby  enhancing  the  repre¬ 
sentativeness  of  the  aerosol  models  used  in  that 
scheme. 

The  total  sea  surface  reflectance  is  then  computed 
for  the  HyperSAS  or  SeaPRISM  viewing  configura¬ 
tion  and  for  the  coincident  aerosol  optical  properties 
data  derived  from  the  collocated  AERONET  system. 
The  determination  of  the  aerosol  mixture  is  based 
on  the  Level  2.0  AERONET  inversion  data  products 
acquired  within  a  time  window  of  plus  or  minus 
20  min  of  the  acquisition  sequence  considered. 
Finally,  based  on  Eqs.  (13)  and  (16),  the  water¬ 
leaving  radiance  is  retrieved  from  the  measured  total 
radiance  Lt  as  follows: 


Lw(O8,0v,4>)  =  L?(08,0„t0)- 


Amrf(0*»  Am  <P)E8(0b) 

K 


(18) 


B.  Normalized  Water-Leaving  Radiance 

The  exact  normalized  water-leaving  radiance,  LWN> 
is  then  obtained  from  Lw(X,0f(p)  as  follows: 


Lwn  =  CtoQ(08,OVf4>,  ra,  IOP,  ws)Cf/Q(08y  ra,  IOP) 


* -w  F* 


(19) 


where  ra  and  IOP  stand  for  the  aerosol  optical  thick¬ 
ness  and  the  inherent  optical  properties,  respec¬ 
tively.  Crq  and  Cf/Q  are  introduced  to  remove  the 
dependence  from  the  viewing  geometry  due  to  the 
refraction  of  wavy  sea  surface  and  the  bidirectional 
effects  in  LWf  respectively;  the  exact  formulations 
of  these  terms  can  be  found  in  [4]  for  example.  The 
term  F0  is  the  theoretical  spectral  extraterrestrial 
solar  irradiance  for  the  mean  Sun-Earth  distance 
[571.2?,  is  downwelling  spectral  irradiance  just  above 
the  sea  surface. 

In  order  to  apply  the  correction  scheme  described 
in  the  previous  section,  the  spectral  values  of  EB 
are  needed.  The  value  of  Es  is  directly  measured 
by  the  irradiance  sensor  of  the  HyperSAS  system. 
In  contrast,  the  SeaPRISM  system  does  not  have  an 
irradiance  radiometer,  andS,  is  accordingly  parame¬ 
terized  in  the  official  SeaPRISM  processing  as 
follows: 


E,(0.)  =  Tm(0')  cos  (0,)D2Fo.  (20) 

where  D2  accounts  for  the  variation  of  the  Sun-Earth 
distance  with  the  day  of  the  year  [58];  Ttot  stands  for 
the  total  atmospheric  irradiance  transmittance, 
which  is  the  sum  of  its  direct  and  diffuse  components. 
In  the  NASA  SeaPRISM  processing,  the  latter  is 
computed  at  a  given  wavelength  X  following  the 
empirical  relationship  [59,60] 
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(21) 

where  ty,  ra,  and  T03  are  the  Rayleigh,  aerosol  and 
ozone  optical  thicknesses,  respectively.  The  Ca  coeffi¬ 
cient  accounts  for  aerosol  phase  function  and  absorp¬ 
tion  properties,  and  its  value  is  typically  around  0.14 
with  a  slight  spectral  dependence.  In  the  SeaPRISM 
processing,  ty  is  computed  taking  into  account  the 
atmospheric  pressure  variations,  ra  is  directly  re¬ 
trieved  by  the  sun  photometer  component  of  the 
SeaPRISM  measurements  for  each  spectral  band, 
and  t03  is  computed  from  the  Total  Ozone  Mapping 
Spectrometer  dataset.  Note  that  the  use  of  the  para¬ 
meterization  of  Eq.  (21)  might  induce  small  discre¬ 
pancies  between  the  HyperSAS  and  SeaPRISM 
irradiance  products  for  the  usual  atmospheric 
conditions  of  the  LISCO  site  [9], 

5.  Application  to  LISCO  Above-Water  Measurements 

In  this  section,  the  proposed  algorithm  for  correction 
of  the  sky  light  reflection  is  evaluated  on  the  basis  of 
the  1.5-year  LISCO  dataset  of  multi-  and  hyperspec- 
tral  measurements.  Specifically,  the  consistency  of 
the  water-leaving  radiance  retrieved  from  the  collo¬ 
cated  instrumentation  is  analyzed  and  the  implica¬ 
tions  of  considering  polarization  are  discussed. 

A.  SeaPRISM  Data 

The  correction  scheme,  described  in  the  previous  sec¬ 
tion,  was  applied  to  the  SeaPRISM  Level  2.0  data 
over  the  1.5-year  period  of  LISCO  acquisitions,  which 
have  been  reprocessed  by  the  AERONET  program  for 
Level  2.0  data  production  [4],  Those  data  were  con¬ 
trolled  to  exclude  any  data  prohibitively  contaminated 
by  sun  glint,  foam,  clouds,  or  any  other  undesirable  en¬ 
vironmental  effects.  Thus,  the  new  correction  scheme 
and  the  standard  one  can  be  meaningfully  compared  to 
infer  how  well  those  schemes  are  handling  reflection  of 
the  diffuse  sky  light  of  the  rough  sea  surface. 

The  LISCO  data  were  processed  with  the  proposed 
algorithm  following  two  modalities:  first  without 
considering  polarization,  that  is  to  say,  the  correction 
was  performed  based  on  scalar  radiative  transfer 
computations.  With  this  modality,  the  water-leaving 
radiances  are  retrieved  for  physical  assumptions  si¬ 
milar  to  those  made  within  the  NASA  processing 
scheme.  The  second  modality  makes  use  of  vector  ra¬ 
diative  transfer  computations,  enabling  it  to  fully 
consider  the  polarization  state  of  the  signal  during 
the  sky-light  correction  step.  The  respective  time 
series  of  exact  normalized  water-leaving  radiances 
retrieved  by  the  standard  and  the  new  correction 
scheme  with  and  without  polarization  consideration 
are  shown  in  Fig.  5.  In  general,  for  all  the  spectral 
bands  except  that  centered  on  413  nm,  the  seasonal 
variations  observed  on  the  three  datasets  are  consis¬ 
tent  over  the  whole  1.5-year  period.  Especially,  the 
distinctive  peak  of  radiance  occurring  after  March 
17,  2010,  resulting  from  an  increase  of  sediment 


concentration  following  a  significant  storm  event 
with  higher  riverine  input  and  water  body  mixing, 
can  be  observed  in  Fig.  5.  Moreover,  the  NASA 
processing  and  the  proposed  algorithm  yield  very 
similar  water-leaving  radiance  retrievals  when  po¬ 
larization  is  ignored.  As  a  result,  it  can  be  concluded 
that  the  seasonal  changes  are  captured  well  by  the 
two  correction  schemes  over  the  1.5-year  datasets 
of  LISCO  acquisitions. 

At  the  longer  wavelengths,  namely  between  491 
and  668  nm,  the  two  modalities  (with  and  without 
polarization)  yield  virtually  similar  results  with  the 
exception  of  the  winter  periods,  where  slight  differ¬ 
ences  can  be  seen  (green  and  red  dots  are  not  super¬ 
imposed  in  Fig.  5).  At  this  period  of  the  year,  the 
Sun  remains  low  on  the  horizon  within  daytime.  The 
amount  of  diffuse  sky  light  is  therefore  high  in 
comparison  to  direct  sunlight.  In  addition,  it  can  be 
remembered  that  the  nondiagonal  terms  of  the 
Fresnel  matrix  (i.e.,  i?*2  and  #13)  exhibit,  for  large 
solar  angles,  values  practically  as  large  as  those  of 
the  R\\  term  (see  Fig.  2),  which  is  usually  used  for 
sky  light  reflection  calculations  in  the  standard  cor¬ 
rection  scheme.  As  a  result,  the  discrepancies  between 
the  standard  and  the  new  polarization-based  correc¬ 
tion  scheme  fire  enhanced  for  large  solar  zenith  an¬ 
gles,  and,  therefor,  the  differences  between  the 
retrieved  water-leaving  radiances  are  more  pro¬ 
nounced  during  winter  time. 

At  the  shorter  wavelengths,  the  optical  thickness 
of  the  air  molecules  is  much  higher  than  at  longer 
wavelengths.  This  optical  behavior  of  the  atmo¬ 
sphere  makes  the  diffuse  skylight  more  intense  and 
polarized  toward  the  shorter  wavelengths.  It  can  be 
readily  seen  in  Fig.  5(a)  that  the  higher  degree  of 
polarization  occurring  at  413  nm,  for  instance,  in¬ 
duces  strong  discrepancies  between  retrievals  when 
polarization  is  taken  into  account  or  not.  Moreover, 
these  discrepancies  are  also  enhanced  in  winter,  as 
discussed  before,  producing  very  large  differences 
between  the  two  correction  schemes. 

This  qualitative  view  offered  by  the  times  series 
enables  us  to  conclude  that  the  polarization  of  diffuse 
sky  light  significantly  impacts  the  water-leaving  re¬ 
trieval  from  above-water  radiometric  measurements 
at  the  shortest  wavelengths.  On  the  other  hand,  a 
more  quantitative  view  on  polarization  impacts  can 
be  reached  by  matchup  comparison  of  the  coincident 
products  obtained  through  application  of  the  two 
different  correction  schemes.  In  order  to  analyze 
these  comparisons  the  following  statistical  indicators 
are  used.  The  unbiased  relative  percent  difference 
(URPD)  and  the  absolute  unbiased  relative  percent 
difference  (|URPD|)  are  defined  similarly  to  [7,38]: 


1  N 

URPD  =  200  x-iY" 

1=1 
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Fig.  5.  (Color  online1*  Time  series  of  the  normalized  water-leaving  radiances  (in  mW  cm'2  sr~*  ^im'1 )  obtained  from  the  LISCO  SeaPRISM 
Level  2.0  data  by  the  NASA  algorithm  (black  dots)  and  the  radiative-transfer-based  algorithm  without  (red  dots)  and  with  polarization 
(green  dots)  at  (a)  413,  (b)  442,  (c)  491,  (d)  551,  and  (e)  668  nm. 


where  x  stands  for  the  standard  products  and  y  for 
the  new  product  to  be  evaluated.  N  is  the  number 
of  matchup  points.  The  root-mean-square  differences 
(RMSD)  complete  those  statistical  indicators.  Finally, 
a  least  squares  fit  is  also  adjusted  within  the 
matchup  points,  with  the  associated  coefficient  of 
determination,  R2f  and  the  equation  of  the  regression 
line.  It  should  be  noted  that  no  statistical  filtering 
has  been  applied  to  the  following  comparisons  to 
remove  outliers. 

The  matchup  comparison  of  the  normalized  water¬ 
leaving  radiances  obtained  from  Level  2.0  SeaPRISM 
data  by  the  standard  and  new  correction  schemes 
when  polarization  is  ignored  is  shown  in  Fig.  6(a) 
for  the  spectral  bands  of  interest.  To  complete  this 
spectral  average  view,  a  summary  of  the  statistical  in¬ 
dicators  calculated  for  each  spectral  band  is  provided 
in  Table  1.  In  this  comparison,  the  two  datasets  were 
built  while  neglecting  polarization  impacts  in  the  re¬ 
spective  data  processing.  In  this  sense,  the  compari¬ 
son  of  Fig.  6(a)  enables  us  to  validate  the  proposed 
algorithm  based  on  radiative  transfer  computation 
and  aerosol  knowledge  because  both  procedures  rely 
on  the  same  physical  assumptions.  On  the  spectral 
average,  the  coefficient  of  determination  R2  is 


almost  1,  showing  a  perfect  correlation  with  a  regres¬ 
sion  line  very  close  to  the  1 : 1  line.  This  strong  corre¬ 
lation  is  observed  for  every  wavelength  as  shown  by 
the  spectral  values  ofi?2  >  0.95  in  Table  1.  The  overall 
dispersion  is  reasonably  low  with  |URPD|  <  5%  and 
an  RMSD  of  0.021  mW  cm**2  sr~*  pm*”1.  This  excellent 
agreement  between  the  two  datasets  demonstrates 
the  correctness  and  reliability  of  the  proposed 
radiative-transfer-based  approach  for  correcting  the 
reflected  sea  surface  signal. 

Because  the  radiative-transfer-based  approach  has 
been  validated  over  the  whole  time  series,  the  pro¬ 
posed  algorithm  can  be  used  to  analyze  the  impact 
of  polarization  on  the  water-leaving  radiance  retrie¬ 
val.  Figure  6(b)  shows  the  comparison  of  the  exact 
water-leaving  radiances  obtained  with  the  proposed 
algorithm  when  polarization  is  neglected  or  fully  con¬ 
sidered.  On  the  spectral  average,  the  coefficient  of 
determination  R2  is  0.99,  showing  a  very  strong  cor¬ 
relation  between  retrievals.  However,  this  parameter 
R2  decreases  rapidly  at  the  shorter  wavelengths. 
At  551  nm,  R2  is  equal  to  0.99  and  decreases  down 
to  0.81  at  413  nm.  At  the  same  time,  the  |URPD|  va¬ 
lues  increases  from  1.7%  to  15.3%  between  551  and 
413  nm  for  a  spectral  average  of  6.5%.  Therefore,  it 
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Fig.  6.  (Color  online)  Comparison  of  the  exact  normalized  water¬ 
leaving  radiances  (in  mW  cm-2  sr-1  pm’1)  derived  from  SeaPRISM 
measurements  through  (a)  the  standard  data  processing  and  the 
radiative-transfer-based  correction  scheme  ignoring  polarization 
and  (b)  the  radiative-transfer-based  correction  modality  ignoring 
polarization  and  that  fully  considering  the  polarization  impacts. 

can  be  concluded  that  the  impact  of  polarization  on 
sky-light  reflection  is  conspicuous  at  the  shorter 
wavelengths  and  is  markedly  less  prohibitive  in  the 
green  and  red  part  of  the  spectrum.  On  the  other 
hand,  it  is  worth  noticing  that  the  URPD,  jURPD|, 
and  RMSD  values  are  in  general  larger  than  those 


observed  in  Fig.  6(a)  and  are  quite  comparable  to 
those  retrieved  in  the  matchup  comparison  between 
above-  and  in-water  measurements  [10,15,38],  prov¬ 
ing  that  the  polarization  effects  related  to  sky-light 
reflection  is  of  major  importance  when  dealing  with 
water-leaving  radiance  retrieval  from  above  water. 

B.  HyperSAS  Data 

In  the  same  manner  as  for  SeaPRISM  data,  the  stan¬ 
dard  and  new  sea  surface-correction  schemes  were 
applied  to  the  HyperSAS  measurements  to  retrieve 
the  exact  normalized  water-leaving  radiance,  Lwat* 
However,  the  setup  of  the  HyperSAS  system  is 
noticeably  different  than  that  of  SeaPRISM  as  the 
HyperSAS  radiometers  are  fixed  on  the  platform, 
and  their  relative  azimuth  with  the  Sun  is  perpe¬ 
tually  changing  over  the  day.  In  order  to  illustrate 
the  efficiency  of  the  respective  standard  and  new  cor¬ 
rection  schemes,  the  spectra  of  LWN  retrieved  within 
a  day  are  displayed  in  Fig.  7.  On  August  19,  2010, 
a  field  cruise  was  conducted  in  the  vicinity  of  the 
LISCO  platform,  enabling  us  to  check  the  clear  sky 
and  low  wind  conditions  of  that  day,  which  was 
chosen  for  the  illustration  of  Fig.  7. 

The  Lwn  spectra  estimated  by  the  standard  correc¬ 
tion  scheme  are  plotted  in  Fig.  7(a)  for  the  different 
relative  azimuths  for  which  the  HyperSAS  measure¬ 
ments  were  performed.  It  can  be  readily  seen  in  this 
figure  that  the  azimuth  variation  induces  slight 
variation  of  the  retrieved  Lyf^  provided  that  the 
azimuth  is  larger  than  70°  for  wavelengths  longer 
than  440  nm,  approximately.  It  should  be  noted  that 
the  quality  check  procedure  of  LISCO  data,  recently 
defined  on  the  basis  of  the  standard  correction 
scheme  [9],  recommended  utilization  of  HyperSAS 
LyjN  only  if  measurements  were  performed  for  azi¬ 
muths  larger  than  70°.  Nevertheless,  the  spectrum 
retrieved  from  measurements  performed  with  an 
azimuth  of  73°  shows  significant  departure  from 
the  spectra  retrieved  for  larger  azimuth,  whereas 
it  would  probably  be  expected  to  measure  a  virtually 
unchanged  spectrum  within  the  same  day  of  interest. 
Furthermore,  when  the  azimuth  is  59°,  the  standard 
correction  fails  to  retrieve  physically  meaningful 
water-leaving  radiance  as  demonstrated  by  the  nega¬ 
tive  values  obtained  for  wavelengths  smaller  than 


Table  1.  Summary  of  the  Statistical  Indicators  for  the  Comparison  of  the  Exact  Normalized  Water-Leaving  Radiances  of  Fig.  6 


Data 

Wavelength  (nm) 

All  Wavelengths 

413 

442 

491 

551 

668 

Figure  6(a) 

R2 

0.951 

0.987 

0.997 

0.999 

0.996 

0.998 

Regression  line 

0.96x  4-  0.01 

0.99*  -  0.01 

1.00* -0.02 

1.01*  -  0.02 

1.01* -0.01 

0.99* -0.01 

URPD 

-0.8 

-5.4 

-2.7 

-1.8 

-3.6 

-2.9 

|URPD| 

7.6 

6.8 

3.2 

1.9 

3.8 

4.7 

RMSD 

2.37  x  10-2 

2.60  x  10-2 

2.16  x  10-2 

1.83  x  lO'2 

1.04  x  10“2 

2.07  x  10~2 

Figure  6(b) 

R2 

0.811 

0.943 

0.989 

0.996 

0.991 

0.991 

Regression  line 

0.95*  +  0.03 

1.02*  -  0.01 

1.02*  +  0.00 

1.01* +  0.00 

1.01*  +  0.00 

1.00* +  0.01 

URPD 

6.3 

3.1 

1.0 

0.5 

0.9 

2.4 

|URPD| 

15.3 

9.1 

3.6 

1.7 

2.8 

6.5 

RMSD 

5.11  x  10-2 

4.30  x  10-2 

3.05  x  lO'2 

2.06  x  10’2 

0.88  xlO"2 

3.43  x  lO’2 
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Fig.  7.  (Color  online)  Exact  watei^leaving  radiance  spectra  (in 
mW  cm'2  sr'1  pnT1)  retrieved  through  (a)  the  standard  and 
(b)  the  new  polarization-based  correction  scheme  applied  to  the 
HyperSAS  data  of  August  19,  2010.  The  colors  indicate  the  differ¬ 
ent  relative  azimuths  <f>  for  which  measurements  were  performed. 

450  nm  and  longer  than  740  nm.  Therefore,  the  LWx 
values  retrieved  by  the  standard  correction  scheme 
are  characterized  by  significant  azimuth  dependence 
even  though  the  measurement  conditions  were  satis¬ 
factory  (e.g.,  low  wind  speed,  clear  sky). 

In  contrast,  application  of  the  new  correction 
scheme  on  the  same  dataset  yields  very  similar  spec¬ 
tra  over  same  the  day  regardless  of  the  azimuth  of 
the  acquisitions;  see  Fig.  7(b).  In  addition,  the  slight 
variations  observed  for  the  intermediate  wave¬ 
lengths  (e.g.,  490-630  nm)  are  most  likely  connected 
to  biological  changes  induced  by  irradiance  variation 
over  the  day.  Regarding  the  infrared  part  of  the  spec¬ 
trum  where  the  optical  impact  of  phytoplankton  is 
negligible,  the  retrieved  L\yn  are  markedly  constant 
for  all  the  azimuths  considered.  The  same  observa¬ 
tion  is  true  for  the  shorter  wavelengths.  Based  on  the 
remarkable  consistency  of  those  HyperSAS  spectra 
within  the  same  day,  it  can  be  safely  concluded  that 
the  new  correction  scheme  yields  meaningful  LWN 
estimations,  thereby  validating  the  polarization- 
based  approach. 

As  for  SeaPRISM  data,  matchup  comparisons 
between  Lwn  retrieved  by  the  standard  and  new  cor¬ 
rection  schemes  applied  to  HyperSAS  data  were 
performed  over  the  1.5-year  dataset.  Figure  8(a) 
shows  the  comparison  between  the  standard  data 
processing  and  the  radiative-transfer-based  correc¬ 
tion  scheme  ignoring  polarization,  whereas  Fig.  8(b) 
shows  the  comparison  between  the  radiative-transfer- 
based  correction  modality  ignoring  polarization  and 


that  fully  considering  the  polarization  impacts. 
The  statistical  indicators  are  summarized  in  Table  2. 
In  those  comparisons,  the  azimuth  range  was  re¬ 
stricted  to  70°-180°  because  the  HyperSAS  retrie¬ 
vals  are  quality  assured  outside  this  range  for  the 
standard  correction  scheme. 

The  comparison  of  the  normalized  water-leaving 
radiances  obtained  from  HyperSAS  data  by  the  stan¬ 
dard  and  new  correction  schemes  when  polarization 
is  ignored  is  shown  in  Fig.  8(a).  In  this  comparison, 
the  two  datasets  were  built  neglecting  polarization 
impacts  in  the  respective  data  processing.  This 
comparison  exhibits  very  similar  results  to  those  of 
Fig.  6(a)  based  on  the  SeaPRISM  datasets.  On  spec¬ 
tral  average,  the  coefficient  of  determination  R 2  is 
almost  1,  showing  a  perfect  correlation  with  a  regres¬ 
sion  line  very  close  to  the  1 1 1  line.  This  strong  corre¬ 
lation  is  observed  for  every  wavelength  as  shown  by 
the  spectral  values  of R2  >  0.97  in  Table  2.  The  overall 
dispersion  is  reasonably  low  with  |URPD|  <  5%  and 
an  RMSD  of  0.026  mWcm^sr"1  pm"1.  Once  again, 
this  excellent  agreement  between  the  two  datasets 


Lws  ( without  polarization) 

Fig.  8.  (Color  online)  Comparison  of  the  exact  normalized 
water-leaving  radiances  (in  mW  cm“2sr‘1  pm-1)  derived  from 
HYPERSAS  measurements  through  the  standard  data  processing 
and  the  polarization-based  correction  scheme.  The  azimuth  range 
is  between  70°  and  180°. 
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Table  2.  Summary  of  the  Statistical  Indicators  for  the  Comparison  of  the  Exact  Normalized  Water-Leaving  Radiances  of  Fig.  8 


Data 

Wavelength  (nm) 

All  Wavelengths 

413 

442 

491 

551 

668 

Figure  8(a) 

R2 

0.971 

0.990 

0.998 

0.999 

0.996 

0.998 

Regression  line 

0.99* -0.03 

0.99.*  -  0.02 

l.OOr-O.02 

1.0 lx  -  0.01 

l.Olx -0.01 

l.OOx  -  0.02 

URPD 

-9.2 

-6.5 

-2.8 

-1.2 

-1.8 

-4.3 

|URPD| 

9.3 

6.6 

2.8 

1.4 

2.4 

4.5 

RMSD 

3.53  x  lO'2 

3.51  x  lO'2 

2.58  x  lO"2 

1.66  x  10~2 

0.85  x  lO'2 

2.64  x  10’2 

Figure  8(b) 

R2 

0.580 

0.865 

0.976 

0990 

0.985 

0.977 

Regression  line 

0.81x  +  0.07 

1.02x-0.01 

1.04x  -  0.03 

l.Olx  -  0.01 

1.02x  -  0.01 

l.OOtr  +  O.OO 

URPD 

1.2 

-0.5 

-0.5 

0.0 

-0.1 

0.0 

|URPD| 

20.6 

12.8 

5.5 

2.8 

4.0 

9.2 

RMSD 

8.53  x  lO'2 

7.56  x  lO'2 

5.45  x  10~2 

3.49  x  lO"2 

1.45  xlO-2 

5.90  x  10-2 

corroborates  the  correctness  and  reliability  of  the 
proposed  radiative-transfer-based  approach  already 
demonstrated  based  on  the  SeaPRISM  datasets. 

Conversely,  the  comparison  between  the  datasets 
generated  in  neglecting  or  accounting  for  polariza¬ 
tion  exhibits  more  pronounced  discrepancies  in 
Fig.  8(b).  On  the  spectral  average,  the  coefficient  of 
determination  of  0.98  shows  a  good  correlation,  but 
the  averaged  |URPDj  value  of  about  9%  and  an  RMSD 
equal  to  0.059  mW  cm”2  sr”1  pm'1  indicate  significant 
dispersion  induced  by  the  polarization  effects.  More¬ 
over,  a  strong  spectral  dependence  is  also  observed 
with  increasing  |URPD|  values  toward  the  shorter 
wavelengths,  namely  2.8%  at  551  nm  and  20%  at 
413  nm,  accompanied  by  a  decrease  of  R 2  from  0.99 
to  0.58.  It  is  worth  noticing  that  all  those  observed  dis¬ 
crepancies  generated  by  the  polarization  state  of  the 
reflected  sky  light  are  generally  more  pronounced  for 
the  HyperSAS  datasets  than  for  the  SeaPRISM  ones. 
This  can  be  explained  by  the  different  geometric 
setups  used  by  the  two  systems:  SeaPRISM  acquisi¬ 
tions  are  performed  for  a  fixed  relative  azimuth  angle 
of  90°  whereas  HyperSAS  data  are  acquired  for  var¬ 
ious  relative  azimuth  angles,  here  comprising  those 
between  70°  and  180°.  As  discussed  in  Section  2.C 
(see  Fig.  4),  the  polarization  effects  on  the  sea  surface 
signal  are  highly  variable  with  azimuth.  In  addition, 
acquisitions  at  small  relative  azimuths  are  more 
impacted  by  the  sky-light  reflection  because  the  dif¬ 
fuse  sky  is  generally  increasing  when  the  sensor  is 
pointing  closer  to  the  solar  plane.  In  conclusion,  the 
polarization  effects  should  be  accounted  for  when  cor¬ 
recting  above-water  measurements  for  the  sea  surface 
signal  (impact  larger  than  9%  on  the  water-leaving  ra¬ 
diance  retrievals)  and  particularly  when  dealing  with 
a  large  range  of  azimuth  angles  as  is  the  case  for  the 
HyperSAS  system  at  the  LISCO  site  or  other  radio¬ 
meters  commonly  in  use  in  oceanographic  field  cruises. 

The  implication  of  considering  polarization  for 
sea  surface-reflectance  estimation  can  then  be  eval¬ 
uated  on  the  basis  of  comparison  of  the  collocated 
and  coincident  SeaPRISM  and  HyperSAS  datasets. 
In  this  case,  comparisons  between  SeaPRISM  and 
HyperSAS  LWN  retrievals  were  carried  out  for  the 
five  main  SeaPRISM  spectral  bands,  namely  413, 
442,  491,  551,  and  668  nm.  The  hyperspectral 


HyperSAS  data  were  integrated  with  the  sensor  re¬ 
lative  spectral  response  function  of  each  SeaPRISM 
band  in  order  to  produce  equivalent  data  for  both 
systems.  The  data  involved  were  restricted  to 
SeaPRISM  measurement  sequences  taken  within 
±10  min  of  HyperSAS  sequence  intervals.  Because 
the  sequence  interval  is  30  min  for  both  systems, 
the  comparisons  were  always  achieved  between  the 
same  single  sequences. 

In  order  to  properly  examine  the  consistency  of 
the  SeaPRISM  and  HyperSAS  datasets,  the  two 
datasets  were  normalized  by  the  same  downwelling 
irradiance  values  measured  by  the  HyperSAS  irradi- 
ance  sensor.  In  the  comparisons,  the  azimuth  range 
of  the  HyperSAS  acquisitions  was  first  considered 
between  70°  and  180°,  corresponding  to  the  good 
data  quality  of  the  HyperSAS  standard  processing. 
The  potential  improvements  of  the  new  correction 
scheme  for  supplementary  azimuth  range,  namely 
azimuth  angles  down  to  30°,  were  then  investigated. 
Results  of  these  comparisons  are  displayed  in  Fig.  9 
for  different  modalities  of  the  water-leaving  radiance 
retrieval  algorithm,  and  statistical  indicators  are 
summarized  in  Table  3. 

In  Fig.  9(a),  the  LWN  retrieved  by  the  standard  cor¬ 
rection  scheme  applied  to  SeaPRISM  and  HyperSAS 
data  are  compared  for  HyperSAS  azimuth  angles 
occurring  between  70°  and  180°.  Notwithstanding 
the  strong  correlation  CR2  =  0.98),  the  comparison 
shows  significant  dispersion  with  |URPD|  =  13.5% 
and  RMSD  =  0.078  mW  cm”2  sr”1  pm”1  on  spectral 
average.  At  413  nm,  where  polarization  impact  is  the 
strongest,  the  respective  datasets  are  less  correlated 
CR2  =  0.81)  and  exhibit  important  dispersions 
between  each  other  (|URPD|  =  24.6%;  RMSD  = 
0.092  mW  cm”2  sr”1  pm”1).  In  contrast,  when  compar¬ 
isons  are  performed  based  on  the  new-correction- 
scheme  retrievals  (Fig.  9(b)),  the  overall  correlation 
is  stronger  with  R2  =  0.98  for  a  |URPD|  value  around 
9.9%.  This  is  particularly  true  at  413  nm  where 
R2  =  0.91  for  a  |URPD|  value  of  12%  and  a  low 
RMSD  of 0.046  mW  cm"2  sr”1  pm”1.  These  values  are 
in  the  same  order  of  magnitude  as  other  studies 
examining  comparison  of  in-water  and  above-water 
datasets;  see  [10]  for  example.  Thus  the  new  correc¬ 
tion  scheme  enables  us  to  reconcile  the  two  datasets 
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Fig.  9.  (Color  online)  Comparison  of  the  exact  normalized  water-leaving  radiances  (in  mW  cm-2  sir1  pm-1 )  derived  from  HYPERS  AS  and 
SeaPRISM  measurements  through  (a),  (c)  the  standard  data  processing  and  (b),  (d)  the  radiative-transfer-based  correction  scheme.  The 
azimuth  range  of  HyperSAS  data  acquisition  is  comprised  (a),  (b)  between  70°  and  180°  and  (c),  (d)  between  30°  and  70°. 


Table  3.  Summary  of  the  Statistical  Indicators  for  the  Comparison  of  the  Exact  Normalized  Water-Leaving  Radiances  of  Fig.  9 


Data 

Wavelength  (nm) 

All  Wavelengths 

413 

442 

491 

551 

668 

Figure  9(a) 

R2 

0.812 

0.921 

0.966 

0.954 

0.953 

0.976 

Regression  line 

0.98x  +  0.08 

0.94x  +  0.10 

0.92x  -1-  0.10 

0.89*  -1-0.11 

0.95r  +  0.04 

0.92x  -1-  0.09 

URPD 

23.8 

14.5 

6.9 

0.4 

11.1 

11.2 

|URPD| 

24.6 

15.8 

8.4 

6.9 

11.6 

13.5 

RMSD 

9.21  x  10~2 

8.43  x  10-2 

7.70  x  10‘2 

8.74  x  10'2 

3.97  x  10"2 

7.84  x  lO'2 

Figure  9(b) 

R2 

0.910 

0.944 

0.970 

0.954 

0.941 

0.980 

Regression  line 

0.92x  +  0.04 

0.9Lc  -1-  0.07 

0.92x  +  0.08 

0.88x  +  0.11 

0.94x  +  0.04 

0.93x  +  0.06 

URPD 

8.0 

8.2 

3.5 

-0.3 

11.4 

6.2 

|URPD| 

12.0 

10.9 

7.1 

7.0 

12.2 

9.9 

RMSD 

4.61  x  10"2 

6.16  x  10 ~2 

7.05  x  lO"2 

8.15  x  10"2 

4.25  x  lO-2 

6.49  x  lO"2 

Figure  9(c) 

R2 

0.086 

0.335 

0.755 

0.778 

0.403 

0.695 

Regression  line 

0.81jc  -  0.08 

1.05*  -  0.14 

1.04jc  -  0.16 

0.93x  +  0.10 

1.07x  +  0.06 

0.94x  -  0.08 

URPD 

-262.7 

-13.1 

-36.3 

-23.5 

20.9 

-62.8 

|URPD| 

305.6 

69.9 

39.1 

23.7 

50.3 

97.7 

RMSD 

4.00  x  10*1 

3.64  x  10-1 

2.73  x  10”1 

2.96  x  10' 1 

1.87  x  10"1 

3.13  x  wr1 

Figure  9(d) 

R2 

0.905 

0.931 

0.970 

0.970 

0.910 

0.976 

Regression  line 

l.OLc  -  0.01 

0.99*  -1-  0.03 

0.95x  -|-  0.02 

0.89tx  -1-  0.08 

0.97jc  +  0.04 

0.90x  -h  0.06 

URPD 

-2.7 

2.8 

-3.5 

-7.2 

10.9 

0.1 

|URPD| 

12.1 

11.3 

8.6 

10.1 

12.6 

11.0 

RMSD 

4.59  x  10-2 

6.64  x  IQ"2 

6.55  x  lO"2 

1.27  x  10“l 

5.45  x  10-2 

7.96  x  lO'2 
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over  the  whole  spectral  range.  It  can  therefore  be 
safely  concluded  that  polarization  effects  on  diffuse 
sky-light  reflection  on  the  sea  surface  must  be  consid¬ 
ered  for  accurate  derivation  of  the  water-leaving 
radiance  from  above  water. 

In  Figs.  9(c)  and  9(d),  the  LWN  retrieved  by  the  stan¬ 
dard  and  new  correction  schemes  are  compared  for 
HyperSAS  azimuth  angles  comprising  between  30° 
and  70°.  It  was  already  noted  elsewhere  that  azi¬ 
muths  smaller  than  70°  might  lead  to  very  low  data 
quality  when  the  standard  processing  is  applied  (9). 
This  behavior  can  be  readily  observed  in  Fig.  9(c), 
where  correlation  between  the  SeaPRISM  data  (ac¬ 
quired  for  a  90°  azimuth  angle)  and  the  HyperSAS 
data  (acquired  for  azimuths  ranging  from  30°  to  70°) 
is  markedly  low  ( R 2  =  0.69),  associated  with  a  prohi¬ 
bitive  dispersion  (jURPD)  ~  98%)  due  to  inaccurate 
correction  for  sea  surface  signal  in  the  HyperSAS 
data  processing.  In  contrast,  the  polarization-based 
processing  applied  to  the  same  raw  data  yields 
successful  retrievals  for  the  same  azimuth  range. 
Figure  9(d)  shows  strong  correlation  (R2  =  0.98)  asso¬ 
ciated  with  satisfactory  dispersion  (jURPDj  <  11%), 
similarly  to  the  results  obtained  for  the  previous 
azimuth  range  of  (70°,180°1.  This  provides  strong  evi¬ 
dence  that  one  of  the  main  limiting  factors  for  water¬ 
leaving  derivation  for  small  azimuth  angles  arises 
from  ignoring  the  polarization  component  of  the  sea 
surface  signal  in  above-water  data  processing.  There¬ 
fore,  it  is  advisable  to  account  for  polarization  in  data 
processing  of  above-water  measurements  when  the 
azimuth  range  is  restricted  to  small  angles  as  can 
happen  in  field  experiments  (e.g.,  due  to  ship  shadow 
avoidance).  To  conclude,  the  spectral  consistency  of 
the  HyperSAS  data  over  the  day  and  the  reconcilia¬ 
tion  of  HyperSAS  and  SeaPRISM  datasets,  thanks 
to  the  polarization-based  correction  scheme,  demon¬ 
strates  the  utility  and  validity  of  considering  the 
whole  physical  phenomenon  (intensity  plus  polariza¬ 
tion  of  light)  of  sky-light  reflection  on  a  rough  sea 
surface. 

6.  Conclusion  and  Perspectives 
The  water-leaving  radiation  measured  from  the 
above-water  level  is  strongly  affected  by  the  effect 
of  the  sky-light  reflection  on  the  ruffled  sea  surface. 
This  contamination  might  be  sometimes  prohibitive 
for  accurate  measurement  of  the  water-leaving  radi¬ 
ance,  Fortunately,  retrievals  of  the  water-leaving  ra¬ 
diance  can  be  achieved  in  most  cases  provided  that 
sun  and  sky  glint  have  been  accurately  handled  and 
subtracted  from  the  total  signal.  In  this  study,  the 
estimation  of  the  sea  surface  reflectance  originating 
from  sky-light  reflection  has  been  investigated  on  the 
basis  of  the  full -ocean -surface  Fresnel  matrix,  which 
includes  reflection  of  the  polarization  terms  of  light 
radiance. 

The  sensitivity  of  above-water  measurements  to 
the  polarization  state  of  the  incident  sky  light  has 
been  theoretically  analyzed  based  on  radiative  trans¬ 
fer  computation.  It  has  been  shown  that  polarization 


of  incident  light  on  the  sea  surface  induces  significant 
variations  of  the  upwelling  radiance  field.  In  particu¬ 
lar,  the  reflectance  of  the  sea  surface  is  strongly 
dependent  on  the  polarization  state  of  the  down  wel¬ 
ling  sky  light,  which  is  mostly  driven  by  the  relative 
proportions  of  air  molecules  and  aerosols.  In  addition, 
the  sea  state  must  be  accurately  handled  to  express 
the  polarization  contribution  to  the  surface  reflec¬ 
tance.  It  was  also  shown  that  errors  in  the  surface 
reflectance  estimation  can  reach  more  than  50%  for 
viewing  geometries  generally  used  in  above-water 
radiometry. 

A  newly  developed  correction  scheme  for  sky-light 
reflection  has  been  described  and  evaluated  based 
on  the  data  acquired  over  1.5  years  at  LISCO,  which 
combines  an  AERONET  SeaPRISM  system  with  a 
hyperspectral  HyperSAS  system.  In  this  scheme, 
the  aerosol  properties  retrieved  by  the  standard 
inversion  algorithm  of  the  AERONET  system  are 
used  to  simulate  the  sea  surface  reflectance  by  in¬ 
cluding  polarization  in  the  resolution  of  the  radiative 
transfer  equation.  Thus,  the  surface  reflectance  is 
retrieved  for  each  multi-  or  hyperspectral  measure¬ 
ment  acquired  from  the  LISCO  platform.  The  nor¬ 
malized  water-leaving  radiance  was  then 

derived  after  subtraction  of  the  sea  surface  compo¬ 
nent  of  the  measured  signal. 

The  validity  and  reliability  of  the  proposed 
radiative-transfer-based  approach  for  correcting 
the  sea  surface  reflected  signal  has  been  confirmed 
based  on  comparison  of  the  Level  2  AERONET-OC 
data  with  the  outputs  of  the  proposed  algorithm 
when  polarization  is  ignored  (i.e.,  R 2  >  0.99  and 
|URPD|  <  5%).  On  the  other  hand,  it  has  been  demon¬ 
strated  that  the  Lwx  retrieved  values  are  strongly 
affected  when  polarization  information  is  taken  into 
account.  The  multispectral  SeaPRISM  data  are  chan¬ 
ged  by  6.4%  and  9.2%  for  the  spectral  average  of  the 
SeaPRISM  and  HyperSAS  system,  respectively.  This 
change  is  greater  than  15%  at  413  nm  where  the  sky  is 
especially  bright.  Generally,  a  stronger  impact  of  the 
polarization  effects  is  observed  for  the  hyperspectral 
HyperSAS  data  because  of  the  fact  that  HyperSAS 
measurements  are  performed  for  variable  azimuths 
during  the  course  of  the  day.  Furthermore,  the  quality 
of  the  HyperSAS  data  was  greatly  improved,  espe¬ 
cially  for  relative  azimuths  smaller  than  70°,  which 
enabled  us  to  expand  the  time  window  of  satisfactory 
HyperSAS  measurements  within  a  day. 

The  collocated  SeaPRISM  and  HyperSAS  Lw N 
datasets  showed  satisfactory  agreement  with  each 
other  when  polarization  was  accounted  for  by  the 
correction  scheme,  whereas  significant  discrepancies 
occur  when  polarization  was  neglected.  In  particular, 
the  relative  differences  between  the  two  datasets  are 
now  lower  than  10%  on  spectral  average  for  a  wide 
range  of  azimuth  angles.  It  was  therefore  concluded 
that  the  polarization  effects  arising  during  sky-light 
reflection  on  the  sea  surface  significantly  impact  the 
water-leaving  radiance  derivation  from  above-water 
measurements  and,  in  turn,  the  quality  of  the  field 
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data  utilized  for  satellite  data  validation  or  potential 
vicarious  calibration  purposes. 

This  polarization-based  algorithm  showed  satis¬ 
factory  performances  for  water-leaving  radiance  de¬ 
rivation  at  the  LISCO  site*  However,  it  needs  to  be 
evaluated  at  other  coastal  sites  exhibiting  different 
optical  characteristics,  or  even  for  above-water  mea¬ 
surements  performed  on  open  ocean  waters.  In  this 
manner,  the  relative  impact  of  polarization  on  above - 
water  radiometry  might  be  quantified  worldwide.  In 
this  case,  the  impact  of  the  different  types  of  aerosols 
should  be  first  analyzed  to  generalize  the  approach. 
On  the  other  hand,  similar  algorithms  as  that  used 
in  this  study  could  be  applied  for  processing  radio- 
metric  datasets  commonly  acquired  with  HyperSAS- 
like  systems  during  field  cruise  experiments  where 
restriction  in  the  possible  viewing  configurations 
might  be  a  supplementary  difficulty  for  accurate  re¬ 
trievals  (e.g.,  ship  shadow  avoidance).  Another  way 
of  improvement  would  be  to  reconsider  the  standard 
wave  slope  distribution  used  to  model  the  roughness 
of  the  sea  surface.  In  this  direction,  ongoing  efforts 
are  already  undertaken  to  analyze  and  model  wave 
dynamics  at  small  time  and  spatial  scales  [61-64] 
which  plays  an  important  role  on  light  reflection 
when  dealing  with  small  field  of  view  and  short  time 
integration  [65].  Nevertheless,  and  in  every  case, 
polarization  effects  at  the  surface  will  have  to  be 
accounted  for. 
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